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We briefly review the present status of nonextensive statistical mechanics. We focus on (i) the 
central equations of the formalism, (ii) the most recent applications in physics and other sciences, (iii) 
the a priori determination (from microscopic dynamics) of the entropic index q for two important 
classes of physical systems, namely low-dimensional maps (both dissipative and conservative) and 
long-range interacting many-body hamiltonian classical systems. 

I. CENTRAL EQUATIONS OF NONEXTENSIVE STATISTICAL MECHANICS 

Nonextensive statistical mechanics and thermodynamics were introduced in 1988 and further developed in 
1991 j|] and 1998 with the aim of extending the domain of applicability of statistical mechanical procedures to 
systems where Boltzmann-Gibbs (BG) thermal statistics and standard thermodynamics present serious mathematical 
difficulties or just plainly fail. Indeed, a rapidly increasing number of systems are pointed out in the literature for 
which the usual functions appearing in BG statistics appear to be violated. Some of these cases are satisfactorily 
handled within the formalism we are here addressing (see for reviews and |^ for a regularly updated bibliography 
which includes crucial contributions and clarifications that many scientists have given along the years). Let us start 
by just reminding the central equations. 

First of all, the exponential function is generalized into the q-exponential function 

e^^ [l + (l-g)x]T^ (qeU). (1) 

We can trivially verify that this (nonnegative and monotonically increasing) function (i) for g — s- 1 yields ef = 
(Vx), (ii) for g > 1, vanishes as a power-law when x — > — oo and diverges at a; = l/{q — 1), and (iii) for q < 1, has a 
cutoff at X = —1/(1 — (7), below which it is defined to be identically zero. If x — *■ we have e^^^l+x (Vg). 
The inverse function of the g-exponential is the q-logarithm, defined as follows: 

x^^'i — 1 

In X = — (g e 7^) . (2) 

1 - g 

Of course Ini x ~\nx (Vx) . If a; ^ 1 we have lug a; ~ x — 1 (Vg) . 
The nonextensive entropic form we postulate is 

5 fc (^p, = l; <ze7^) , (3) 

where W is the total number of microscopic configurations, whose probabilities are {pi\. Without loss of generality 
we shall from now on assume k ~ 1. We can verify that, for q — )■ 1, this entropy reproduces the usual Boltzmann- 
Gibbs-Shannon one, namely 5*1 = — X^i^i-P* l^iK- The continuous and the quantum expressions of 5*^ are respectively 
given by 

_ I- Sdx[p{x)]'i 

- ^ (4) 



and 
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where p is the matrix density. Unless specifically declared in what follows, we shall be using the form of Eq. (3). It 
is easy to verify that all its generic properties can be straightforwardly adapted to both the continuous and quantum 
cases. 



Sq can be written as 



Sq = (K (6) 
Pi 



where the expectation value ((...)) = "Y^jiLiPii---)- K can also be written as 

Sq= {-\nqPi)q , (7) 

where the unnormalized q-expectation value is defined to be = (•■■)■ Of course — ((••■))• This is 

a good point for defining also the normalized q-expectation value 

(((■■■))), ^ ^^^^'V\ (8) 

T,^=lPi 

which naturally emerges in the formalism. We verify trivially that = = {{■■■)), and also that {{{■■■)))q — 

((••■)>./(l>,- 

If A and B are two independent systems (i.e., p^^^ — P^pf j)), then we have that 

Sq{A + B)= Sq{A) + Sq{B) + (1 - q) S q{A) S q{B) . (9) 

It is from this property that the name nonextensive statistical mechanics was coined. The cases q < 1 and q > 1 
respectively correspond to superextensivity and subextensivity of Sq since in all cases Sq > 0. 
At equiprobability, i.e., pi = 1/W, we obtain 

Sq ^lHqW , (10) 

which is the basis for the microcanonical ensemble. 

For the thermal equilibrium corresponding to the canonical ensemble of a Hamiltonian system, we optimize Sq with 
the constraints ^^iPi = 1 and {{ei))q — Uq, where {ci} are the eigenvalues of the Hamiltonian of the system, and 
Uq is the generalized internal energy. We obtain |^ 

: : — : : ^ — 7^ > (n) 

^1 



9 [l + (g-l)/3(e,-C/,)]— [l + (g-l)/3'e,] — 

where P is the Lagrange parameter, Zq = ^'^'^^ ^''^ P' ^ ^^^^ defined function of /3. For q = 1 we recover 

the celebrated BG weight. When /3 > and the energy increases, the probability decays like a power law for g > 1 
and exhibits a cutoff for g < 1. 

Analogously, if we optimize 5"^ as given by Eq. (4) with the constraints J dx p{x) — 1 and {{x'^))q ^ {a > 0), we 
obtain the ^-generalization of the Gaussian distribution, namely 

Jdyeq^y [l + (g-l)/32;2]<j-i 

where (3 can be straightforward and explicitly related to a. The variance of these distributions is finite if g < 5/3 and 
diverges if 5/3 < g < 3. For q = 2 we have the Lorentzian distribution. For q > 3 the function is not normalizable, 
and therefore is unacceptable as a distribution of probabilities. 

Let us now address typical time dependences. Let us assume that ^(i) is a quantity characterizing an exponential 
behavior and satisfying ^(0) = 1. Such is the typical case for the sensitivity to the initial conditions of a one- 
dimensional chaotic system, where ^(t) = lim^^(Q')^Q[Ax(t)/Ax(0)], where Ax{t) is the discrepancy at time t of two 
trajectories which, at i = 0, started at a;(0) and at x'{0). Another example is a population which relaxes to zero. If 
N{t) is the number of elements, then we can define ^(t) = N{t)/N{0). The quantity £^{t) monotonically increases in 
our first example (sensitivity), whereas it decreases in the second one (relaxation). The basic equation that ^ satisfies 
is generically 
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(13) 



hence £,{t) — e^^*. In our example of the chaotic system, Ai is the Lyapunov exponent. In our population example, 
we have Ai = — 1/ti, where ti is the relaxation time. What happens in the marginal case Ai = 0? Typically we have 

e = A, e% (14) 

hence 

^ = e^«* = [l + (l-g)A,t]T^ = - — ^. (15) 

[l + {q-l){-X,)t] — 

This quantity monotonically increases if Ag > and q < I, and decreases if Xq = —\/Tq < and q > 1. In both cases 
it does so as a power law, instead of exponentially. In the limit t ^ 0, we have ^ ~ 1 + A, t (Vg). 

A more general situation might occur when both Ai and Xq are different from zero. In such case, many phenomena 
will be described by the following differential equation: 

C = Ale + A,e^ (16) 

hence 



Ag ^ {l-q)Xit 

Ai Ai' 



(17) 



If g < 1 and < Ai << A^, ^ increases linearly with t for small times, as t'^-i for intermediate times, and like e^^* for 

large times. If q > 1 and > Ai >> A<j, ^ decreases linearly with t for small times, as l/t~ for intermediate times, 
and like e~'^^'* for large times. 



II. APPLICATIONS IN AND OUT FROM EQUILIBRIUM 

A considerable amount of applications and connections have been advanced in the literature using, in a variety 
of manners, the above formalism. They concern physics, astrophysics, geophysics, chemistry, biology, mathematics, 
economics, linguistics, engineering, medicine, physiology, cognitive psychology, sports and others It seems ap- 
propriate to say that the fact that the range of applications is so wide probably is deeply related to and reflects the 
ubiquity of self-organized criticality fractal structures and, ultimately, power laws in nature. In particular, a 
natural arena for this statistical mechanics appears to be the so called complex systems . 

We shall briefly review here four recent applications, namely fully developed turbulence ||l0|-[l^, hadronic jets pro- 
duced by electron-positron annihilation jl3j , motion of Hydra viridissima [|l^ , and quantitative linguistics ||l^ . 

Fully developed turbulence: 

As early as in 1996 Boghosian made the first application of the present formalism to turbulence |l^. That was for 
plasma. What we shall instead focus on here is fully developed turbulence in normal fluids. Ramos et al advanced in 
1999 ijl^ the possibility of nonextensive statistical mechanics being useful for such systems. The idea was since then 
further developed by Beck ||T^ and by the Arimitsu's [Q, basically simultaneous and independently. They proposed 
theories within which the probability distribution of the velocity differences and related quantities are deduced from 
basic considerations. We present in Fig. 1 the comparison of Beck's theoretical results with recent high precision 
experimental data for Lagrangian and Eulerian turbulences jl^ . In Fig. 2 we show an analogous comparison between 
Arimitsu's theoretical results and recent computer experimental data. 

Hadronic jets: 

High energy frontal collisions of electron and positron annihilate both and produce relativistic hadronic jets. The 
distribution of the tranverse momenta of these jets admits, as advanced by Fermi, Feynman, Hagedorn and others, a 
thermostatistical theoretical approach. Hagedorn's 1965 theory was q-extended by Bediaga et al in 1999. Their results 
[Q, as well as related ones by Beck [Q, compare quite well with the available CERN experimental data, as can be 
seen in Fig. 3. As important as this is the fact that both Bediaga et al and Beck theories recover a crucial feature of 
Hagedorn's scenario, namely that the temperature to be associated with the distributions should not depend on the 
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collision energy. 
Hydra viridissima: 

The motion of cell aggregates of Hydra viridissima in physiological solution was studied by Upadhyaya et al [ p^ . 
The strongly nonmaxwellian distribution of velocities is quite well fitted with q ~ 1.5: See Fig. 4. They also carried 
diffusion measurements and verified that diffusion is anomalous. Under the assumption that it is of the correlated 
type addressed in [|l9|, they obtained once again q ~ 1.5: See Fig. 5. In other words, two different experiments of 
motion were fitted by one and the same value for q. 

Linguistics: 

The frequency p of words used in a text (say a book or a set of books of one or more authors) as a function of 
their rank r roughly follows the so called Zipf law, namely p cx I/r. This law was improved by Mandelbrot in the 
form p oc l/(a + r)'' . This form fits better real data and is sometimes called the Zipf- Mandelbrot law. This law 
precisely is Eq. (15) with the notation (t,^) = {r,p) (or equivalently Eq. (11) where the energy plays the role of the 
word rank), as first argued by Denisov in 1997 However, although quite faithful at low and intermediate word 
ranks, the Zipf-Mandelbrot law fails at high ranks. This point has been addressed recently by Montemurro [^: See 
Fig. 6, where results are shown from well known texts in English, Greek, Italian and Spanish, from authors of very 
different historical periods and literary styles (all reasonably well fitted by using q ~ 1.9 hence ^ ^ l/{q — I) ~ 1.1). 
Montemurro has shown that the fittings are all sensibly improved by using the present Eq. (17) instead of Eq. (15). 
The agreement obtained is illustrated with a large set of books by Dickens as shown in Fig. 7. In other words, at 
large values of the word rank a crossover occurs from g~1.9to(7~l. The reason for this interesting phenomenon is 
unknown. As a plausible hypothesis, we would like to advance that it might be related to the fact that most authors 
possibly use the very rare words in a manner which reflects their relatively poor knowledge of their exact meaning. 
This attitude could make those words to be used slightly uncorrelated with the context within which they are placed. 
It is however clear that this phenomenon is a very subtle one, and its full elucidation would presumably require very 
sophisticated analysis. 



III. A PRIORI DETERMINATION OF THE ENTROPIC INDEX Q 

Physical bounds to the admissible values of q were first discussed as early as 1993 by Plastino and Plastino (for 
self-gravitating systems), and since then by many others. However, the strict aim of the present Section is how is to 
be determined the value of q to be associated with a specific system whose microscopic (or mesoscopic) dynamics is 
exactly known. This is to say, how the knowledge of the rules that provide the (continuous or discrete) time evolution 
of the system can be used in order to determine without ambiguity the appropriate value(s) of q for that system. This 
crucial question must be answered for the present proposal to be a complete theory, in the sense that it is in principle 
able to predict the results to be expected in all types of experiments with well defined systems. This question has 
by now been answered in several important classes of systems. We shall briefly review here two of them, namely 
low-dimensional maps and long-range many-body classical Hamiltonian systems. 



A. Low-dimensional maps 

We shall focus on one- and two-dimensional maps. The one-dimensional maps necessarily are dissipative. The 
two-dimensional ones can be either dissipative or conservative, but we shall primarily address the latter. Indeed, on 
one hand the dissipative two-dimensional maps are dynamically not so different from the one-dimensional ones. On 
the other hand, two-dimensional conservative maps provide us an avant premiere of Hamiltonian systems. 

One- dimensional maps: 
We consider here one-dimensional dissipative maps of the type 

xt+i ^ f{xua,z) (t = 0, 1,2, e [a;™„,a;„aa;]) , (18) 

where a € 7?. is a control parameter such that when it increases for fixed z, it makes the map to become chaotic (we 
note Oc the smallest value of a above which the system can be chaotic; z G 7?. is another control parameter which differs 
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essentially from a in the sense that z controls the universality class of the chaotic attractor emerging at a = adz)] 
the function / is such that chaotic and nonchaotic behaviors are possible for the variable x, depending on the values 
of (a, z). A paradigmatic such map is the so called z-logistic map, defined as 

xt+i = l-a|xt|" = 0,1,2,...; e hi, 1]) , (19) 

with a G (0,2] and z > 1. The critical value ac{z) {chaos threshold or edge of chaos) monotonically increases from 1 
to 2 when z increases from 1 to infinity; ac(2) = 1.401155... . For z = 2 this map is, as well known, isomorphic to 
Xt+i en Xt{l — Xt). The z-logistic maps and several others have already been studied |Q within the nonextensive 
scenario. We briefly review here their main properties. Most of these properties have been found heuristically, and 
no theorems or rigorous results are available. Consequently, we are unable to precisely specify how generic are the 
properties we are going to describe. We know, however, that wide classes of maps do satisfy them. 

Let us first address the sensitivity to the initial conditions. For all values of a for which the Lyapunov exponent Ai 
is nonzero we verify that q = 1, i.e., £^{t) — *, with Ai < for most values of a < Oc, and Ai > for most values 
of a > flc. However, for the infinite number of values of a for which Ai = we verify that q ^ 1- More precisely, 
for values of a such as those for which bifurcations occur between finite cycle attractors of say the z-logistic map, we 
verify the validity of Eq. (15) with q > 1 and Aq < (this has been very recently proved ||2^). For a — adz) we 
verify that ^(i) exhibits a complex behavior which has, nevertheless, a simple upper bound which satisfies Eq. (15) 
with q < 1 (from now on noted qsen{z), where the subindex sen stands for sensitivity) and Xq{z) > (also this has 
been very recently proved [p4[) . For the universality class of the z-logistic map we verify that qsen monotonically 
increases from minus infinity to a value slightly below unity, when z increases from 1 to infinity (gseri(2) — 0.2445...). 
For the z-cercle and other maps we verify similar behaviors. 

Let us now address the attractor in x space for a = ac(z). Its anomalous geometry can be usefully characterized 
by the so called muUifractal function f{a,z) which typically is defined in the interval amin{z) < a < Oimax{z), and 
whose maximal value is the fractal or ffawsdor/f dimension df{z). For the z-logistic map universality class we have 
d/(z) < 1, whereas for the z-circle map universality class we have d/(z) = 1 (Vz). In all the cases we have checked, 
we verify a remarkable scaling law, namely [ p5[ 

^ ^ ^ (Vz)- (20) 



1 ^sen(z) ^mini,^) ^inax{^) 

This relation has purely geometrical quantities at its right hand member, and a dynamical quantity at its left hand 
member. It can also be shown that 

1 . . In a p ., , . , . 

(^-1)^ (V^), (21) 



l-qsen{z) In 6 

where ai? is one of the two well known Feigenbaum constants, and h is the attractor scaling (6 = 2 for period-doubling 
bifurcations; b = 2/(\/5 — 1) for cercle maps). 

Let us now address the entropy production dSq{t)/dt. We first make a partition of the interval [xmin, Xmax] into 
W nonoverlaping little windows. We place (randomly or not) inside one of those W windows a large number N of 
initial conditions, and run the map t times for each of these points. We generically verify that the N points spread 
into the windows, in such a way that we have the set {Ni}{t) (with ^i{t) = N yt ) . With these numbers we 

can define the set of probabilities {pi{t)} where Pi{t) = Ni{t)/N (Vi). We then choose a value for q and calculate 
Sq{t) by using Eq. (3). We then make an average {Sq){t) over a few or many initial windows (see |2|,|2^ for details), 
and finally evaluate numerically limf_»oo hniw^oo liiiiAr^oo(5'<j)(t)/t. We verify a very interesting result [p6| , namely 
that this limit is finite only for q — qsen{z)', it diverges for all q < Qseniz) and vanishes for all q > qsen{z). We shall 
note this limit Kq and constitutes a natural ^-generalization of the Kolmogorov-Sinai entropy. Summarizing, 

Kq = lim lim lim . (22) 

It is easy to verify, whenever Ai > 0, that qsen{z) — 1 and that the Pesin identity holds, i.e., Ki = Ai. A fascinating 
open question constitutes to find, whenever Ai = (more specifically for a = adz)), under what circumstances the 
conjecture Kg^^^ = Xq^^^ could be true, Xq^^^ being the coefficient appearing in Eq. (15) for q = qsen- We would then 
have the generalization of the Pesin identity. 

Let us next address another aspect |^ concerning the edge of chaos ac(z). We spread now, at t = 0, the N points 
uniformly within the entire [xmin,Xmax] interval, i.e., over all the W windows, and follow, as function of time t, the 
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shrinking of the number W{t) of windows which contain at least one point (disappearence of the Lebesgue measure 

on the X-axis); 14^(0) = W. It can be verified that, for the sequence hmvK->cx3 hmjv^oo, we asymptotically have 

1 

W{t) oc l/t'''-<!i''>-^ , where qeq{z) > 1 (the subindex rel stands for relaxation). The entropic index qrei monotonically 
increases when z increases from 1 to infinity; also, within some range it is verified | p7[ | that l/[qrei{z)~l] oc [l — df{z)]'^. 
We shall now advance a recently established p8[ relation between Qrei and qsen ■ 

Let us go back to the procedure when, sd t ~ 0, only one among the W windows is populated. That single window is 
chosen to be that which makes Sq^^^ {t) to achieve the highest value as t increases. To be more precise Sq^^^ (0) = and 
< Sq^^^{oo) < \nq^^^ W. While t increases, there are many windows for which Sq^^^{t) overshoots above Sq^^^{oo). 
We choose the window for which the overshooting is the most pronounced. After Sq^^^ {t) achieves this peak, it relaxes 

slowly towards ^^^^.^(oo). It does so as l/f^i''-^-^^^'- , where qreiiz,W) approaches its limiting value qrei{z, oo) while 
W diverges. The remarkable fact is that qrei{z, oo) = qrei{z) ! More than this, the approach is asymptotically as 
follows: 

qTel{z) - qrel{z, W) OC ^J^^(^) ■ (23) 

This relation is a remarkable connection between the mixing properties, the equilibration (or relaxing) ones, and the 
degree of graining (from coarse to fine graining while W increases). We may also say that in some sense Eq. (23) 
provides a connection between the Boltzmannian and the Gibbsian approaches to statistical mechanics. Indeed, the 
concept of qsen is kind of natural within a typical Boltzmann scenario where individual trajectories in phase space 
are the "protagonists of the game" , whereas qrei is kind of natural within a typical Gibbs scenario where the entire 
phase space is to be in principle occupied. Before taking into Eq. (23) the qsen — qrei = 1 particular case (i.e., the 
BG-like case), some adaptation is obviously needed; as written in Eq. (23), it is valid only for qsen < 1 and qrei > 1- 



Two-dimensional maps: 
We consider here two-dimensional conservative maps of the type 



xt+i = fx{xt,yt;a,z) 

yt+1 = fy{xt,yt;a,z) (24) 

where xt € [xmin, Xmax] and yt G [ymin,ymax] with t ~ 0,1,2,...; the control parameter z characterizes, as for the 
one-dimensional maps we considered above, the universality class; the control parameter a > and we assume that, 
while a increases from zero to its maximum value, the nonnegative Lyapunov exponent Ai monotonically increases 
from zero to its maximum value. Since the map is conservative (i.e., \d{xt+i, yt+i) / d{xt, yt)\ — 1), the other Lyapunov 
exponent is — Ai. A paradigmatic such map is the so called standard map, defined as follows 

yt+i ^yt + sin(27ra;t) {mod 1) 

Xt+1 = yt+1 + Xt ^ yt + sin(27rxt) + Xt {mod 1) (25) 
Ztt 



as well as its z-generalization |29|, defined as follows 

a 

2^ 



yt+1 = + TT- sin(27ra;t)| sin(27ra;t)|^ ^ {modi) 



xt+i = yt+1 + Xt ^ yt + sin(27rxt)| sin(27ra;t)|^ ^ -I- Xt {mod 1) , (26) 
Ztt 

where z G TZ. 

Some (not clearly characterized yet) classes of such maps exhibit for the entropy production dSq{t)/dt a behavior 
which closely follows the crossover behavior associated with Eq. (17). Let us be more precise. We first partition the 
accessible {x,y) phase in W nonoverlaping little areas (for example W little squares), and put a large number N of 
initial conditions inside one of those areas. As before, we follow along time the set of probabilities {pi}, with which 
we calculate Sq{t) for an arbitrarily chosen value of q. We then average over the entire accessible phase space and 
obtain {Sq){t). Finally we numerically approach the quantity Sq(i) = \\m.w^oo ^^'^N^oo{Sq){t). 

For large values of a, we verify pO| that Si(t) asymptotically increases linearly with i, as expected from the 
fact that Ai(a) > 0, in agreement with Pesin identity. However, an interesting phenomenon occurs for increasingly 
small a, hence increasingly small Ai. For small t (say < i << ti{a,z)), Sg(t) is linear with t for q = 0, and 
acquires an infinite slope for any q < 0. For intermediate t (say ti{a,z) « t << t2{a,z)), Sq{t) is linear with t for 
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q = qsen{z) < 1, acquires an infinite slope for q < qsen{z) and acquires a vanishing slope for q > qsen{z). For large t 
(say t >> t2{a, z)), Sq(t) is linear with t for q = I, acquires an infinite slope for q < 1 and acquires a vanishing slope 
for q > 1. The characteristic times ti{a,z) and t2{a,z) respectively correspond to the [q = 0] ^ [q ^ qsen{z)] and 
[q — qseniz)] — > [g = 1] crossovers. The remarkable feature is that, in the limit a — > 0, ti(a, z) remains finite whereas 
t2{a,z) diverges. In other words, for asymptotically small values of a, the time evolution of Sq(t) is, excepting for 
an initial transient, basically characterized by qseniz) < 1. This fact opens the possibility for something similar to 
occur for Hamiltonian classical systems for which the Lyapunov spectrum tends to zero. This is precisely what occurs 
when the size of the system increase in the presence of long-range interactions, as we shall see in the next Subsection. 
Before closing this subsection, let us mention that studies focusing on qrei for conservative maps are in progress. 

B. Long-range many-body classical hamiltonian systems 

From the thermodynamical viewpoint it is interesting to classify the two-body interactions (and analogously, of 
course, the many-body interactions). According to their behavior near the origin, i.e., for r ^ 0, potentials could be 
classified as collapsing and noncollapsing. Collapsing are those which exhibit a minimum at r = 0. This minimum 
can be infinitely deep, i.e., the potential can be singular at r = 0; such is the case of attractive potentials which 
asymptotically behave as — l/r"^ with > (e.g., Newtonian gravitation, hence v — 1). Alternatively, the potential 
at this r = minimum can be finite, as it is the case of those which behave as —a + br~'^ with a > 0, & > and 
!/ < 0. Collapsing potentials, especially those of the singular type, are known to exhibit a variety of thermodynamical 
anomalies. Noncollapsing potentials are those which exhibit a minimum either at a finite distance (e.g., the Lennard- 
Jones one, or the hard spheres model or any other model having a cutoff at a finite distance vq) or at infinity (e.g., 
Coulombian repulsion). Potentials can be also classified according to their behavior at r — > cx). We can divide them 
into short- and ^orif^-range interactions. Short-range interactions are those whose associated force quickly decreases 
with distance, for example potentials which exponentially decrease with distance, or classical potentials of the type 
— 1/r" with a > d, d being the space dimension where the system is defined. For classical systems, thermodynamically 
speaking, short-iSLUge interactions correspond to the potentials which are integrable at infinity [Q, and long-ia,nge 
interactions correspond to those which are not integrable in that limit, such as those increasing like l/r" with a < 
(which belong to the confining class of potentials, i.e., those which make escape impossible), or those like — l/r" with 
< a/d < 1 (which belong to the nonconfining class of potentials, i.e., those which make escape possible). Long- 
range interactions, especially those of the nonconfining type, are also known to induce a variety of thermodynamical 
anomalies. From the present standpoint a particularly complex potential is Newtonian gravitation (corresponding to 
a — I and d — 3). Indeed, it is both singular at the origin, and long-ranged since < a/d = 1/3 < I. 

In this Section we address an important case, namely that of nonsingular attractive long-range two-body interactions 
in a d-dimensional A'^-body classical hamiltonian system with A^ >> 1. Such systems are being actively addressed 
in the literature by many authors (see and references therein). As an illustration of the thermodynamical 

anomalies that long-range interactions produce, we shall focus on the d-dimensional simple hypercubic lattice with 
periodic boundary conditions, each site of which is occupied by a classical planar rotator. All rotators are coupled 
two by two as indicated in the following Hamiltonian: 

^-Ef +E '''°t^''^^ (^.G[0,2.];a>0). (27) 

j=l 

The distance (in crystal units) between any two sites is the shortest one taking into account the periodicity of the 
lattice. For d = 1, it is r,y = 1,2,3,...; for d = 2, it is Vij = l,\/2, 2,...; for d = 3, it is r^j — 1, \/2, \/3, 2, and 
so on for higher dimensions. We have written the potential term in such a way that it vanishes in all cases for the 
fundamental state, i.e., 6i = (V«), where, without loss of generality, we shall consider 6a — for simplicity. Also 
without loss of generality we have considered unit moment of inertia and unit first-neighbor coupling constant. It 
is clear that, excepting for the inertial term, the present model is nothing but the classical XY ferromagnet. The 
cases a = and a ^ oo respectively correspond to the so called HMF model |3^] (all two-body couplings have 
the same strength), and to the first-neighbor model. For BG statistical mechanics to be applicable without further 
considerations, it is necessary that the potential be integrable, i.e., J-^ dr j-d-^^-a ^ This implies a > d. In this 
case, the energy of the system is extensive, i.e., the energy per particle is finite in the thermodynamic limit A^ — s- oo. 
But the situation becomes more delicate for < a/c? < 1, since then that integral diverges. However, strictly speaking, 
the system being finite, the integral that is to be analyzed is not the one already mentioned but the following one 
instead: 
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drr^-W'", (28) 



which, in the N oo hmit, converges for a/d > 1 and diverges otherwise. It is in fact convenient to introduce the 
quantity 

N=l + d dr r^-^r-" = . (29) 

J I I -a/d 

N equals for a = 0, and, for A'' — > oo, diverges hke N^~"'/'^ for < a/d, diverges hke InTV for a/d = 1, and is 
finite for a/d > 1, being unity in the hmit a/c? — > oo . In general, the energy per particle scales with iV; in other 
words, the energy is nonextensivc for < a/d < 1. To make the problem artificially extensive even for a/d < 1, the 
Hamiltonian can be written as follows: 

The rescaling of the potential of this model is more properly taken into account by ^ij"' @1 rather than by TV, 

but since their N ^ oo asymptotic behaviors coincide, we can as well use N as introduced here. The original (Eq. 
(27)) and rescaled (Eq. (30)) versions of this model are completely equivalent (see [Q) and lead to results that can be 
easily transformed from one to the other version. To make easier the comparison of results existing in the literature, 
we shall from now on refer to the rescaled version (30) . 

The a = model (HMF) clearly is d-independent and is paradigmatic of what happens for any a such that 
< a/d < 1. When isolated (microcanonical ensemble) the a/d = model exhibits a second-order phase transition 
at u = U/N — 0.75, where U is its total energy and N oo. This critical value Uc smoothly increases with 
a/d approaching unity. Dynamical and thermodynamical anomalies exist in both ordered and disordered phases, 
respectively for u < Uc and u > Uc- Let us discuss some anomalies for u > Uc, then some for u < Uc, and finally show 
that these anomalies on both sides of Uc are in fact connected. 

The Lyapunov spectrum is made by couples of real quantities that are equal in absolute value and opposite in sign, 
whose sum vanishes in accordance with the Liouville theorem. The sum of the positive values equals the Kolmogorov- 
Sinai entropy, in accordance with the Pesin theorem. If the maximal Lyapunov exponent vanishes, the entire spectrum 
vanishes, and no exponentially quick sensitivity to the initial conditions is possible. 

We address first the case u > Uc- For the Hamiltonian of rotors we are interested in (Eq. (30)), the maximal 
Lyapunov exponent A™*^^ scales, for large N , like 

Xmax ^(u, a, d) 

where l{u,a,d) is some smooth function of its variables, and K{a/d) decreases from 1/3 to zero when a/d increases 
from zero to unity; k remains zero for all values of a/d above unity, consistently with the fact that A"*"^ is positive 
in that region. In other words, above the critical energy, the sensitivity to the initial conditions is exponential for 
a/d > 1, and subexponential (possibly power like) for < a/d < 1. 

We address now the case < u < Uc, and focus especially on the region slightly below the critical value Uc (e.g., 
u — 0.69 for the a = model). For < a/d < 1, at least two (and possibly only two with nonzero measure) 
important basins exist in the space of the initial conditions: one of them contains the Maxwellian distribution of 
velocities, the other one contains the water-bag (as well as the double water-bag) distibution of velocities. When the 
initial conditions belong to the Maxwellian basin, the system relaxes quickly onto the BG equilibrium distribution 
(strictly speaking, we do not have this numerical evidence but a weaker one, namely that the marginal probability of 
one-rotator velocities tends to the Maxwellian one when N — > oo). When the initial conditions belong to the other 
basin, it first relaxes quickly to an anomalous, metastable (quasi-stationary) state, and only later, at a crossover time 
T, starts slowly approaching the BG equilibrimn. The crossover time diverges with N for a = |^^. It has been 
conjectured that it might in general diverge like t iV. It has been recently established |^ that, for d — 1 
and fixed N , t exponentially vanishes with a approaching unity. All these features are consistent with the conjecture 
T ^ N, which might well be true. During the metastable state, the one-particle distribution of velocities is clearly non 
Gaussian, and in fact it seems to approach the distribution of velocities typical of nonextensivc statistical mechanics 
for g > 1. This anomalous behavior reflects on the sensitivity to the initial conditions. The maximal Lyapunov 
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exponent A"*°^ remains during long time, in fact until t '--^ r, at a low value and then starts approaching a finite value. 
This low value scales like The remarkable feature which has been observed |^,^ for the d = 1 model is 

that k' = k/3 for all values of a. The anomalies above and below the critical point become thus intimately related. 

The whole scenario is expected to hold for large classes of models, including the classical n- vector ferromagnetic-like 
two-body coupled inertial rotors {n = 2 being the present one, n = 3 the Heisenberg one, n ^ oo the spherical one, 
etc). For all of them, in the isolated situation, we expect (i) at the disordered phase, that the maximal Lyapunov 
decreases with TV with the exponent K{a/d,n) (it is yet unclear whether this exponent depends on n or not); (ii) 
at the ordered phase, and starting from initial conditions within a finite basin including the water-bag, that a 
metastable state exists with non BG (possibly g-type) distribution associated with a maximal Lyapunov exponent 
which decreases with N with the exponent K{a/d)/3. In these circumstances, for a/d < 1 (nonextensive systems), the 
limAT^oo limt^oo ordering is expected to yield the usual BG equilibrium, whereas the limj^oo limjv^oo ordering yields 
a non BG (meta)equilibrium, possibly of the type predicted by nonextensive statistical mechanics. This interesting 
phenomenon disappears for a/d > 1 (extensive systems); indeed, both orderings lead then to the same equilibrium, 
namely the BG one, as known since long. 

IV. CONCLUSIONS 

We have presented some of the main peculiarities associated with nonextensive systems. Most of the paradigmatic 
behaviors are expected to become (or have been shown to become) power laws instead of the usual exponentials: 

(i) the sensitivity to the initial conditions is given by ^ = eq' * (typically q < 1); 

(ii) the finite entropy production (Kolmogorov-Sinai entropy like) occurs only for Sq (with q < 1, the same as 
above); 

(iii) the relaxation towards quasi-stationary or (metaequilibrium) states, or perhaps from these to the terminal 
equilibrium states, may occur through Cq (typically q > 1); 

(iv) the stationary, (meta)equilibrium distribution for thermodynamically large Hamiltonian systems may be given 
by Pi cx (typically q > 1, possibly the same as just above). 

The two-dimensional conservative maps exhibit, in the vicinity of integrability and at intermediate times, features 
very similar to those observed in one-dimensional dissipative maps at the edge of chaos. The intermediate stage has 
a duration which diverges when the control parameters approach values where the system is close to integrability. 
Isolated classical Hamiltonian systems behave similarly to low-dimensional conservative maps, 1/iV playing a role 
analogous to the distance of the control parameters to their values where integrability starts. 

The scenario which emerges is that sensitivity and entropy production properties are related to one and the 
same value of qsen < 1 (also related to amin and amax of some multifractal function), whereas the relaxation and 
(meta)equilibrium properties are related to (possibly one and the same) value of g^eZ > 1 (also related to the Haus- 
dorff dimension of the same multifractal function). These two sets of properties are quite distinct and generically 
correspond to distinct values of q (namely, qsen and qrei)- It happens that for the usual, extensive, BG systems they 
coincide providing q^en = qrei = 1, which might sometimes be at the basis of some confusion. In all cases, once the 
microscopic dynamics of the systems is known, it is in principle possible to determine a priori both qs^n and qrei (as 
well as the connection among them and with the chosen graining, as illustrated in Eq. (23)). We have here shown 
how this is done for simple systems. This type of calculation of q from first principles has also been illustrated for a 
variety of other systems The quest for such calculations for more complex systems is in progress. 
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FIG. 1. The distributions of difFerences of velocity in two different experiments of fluid turbulence. The solid lines correspond 
to Beck' s theory (it). 

FIG. 2. The distributions of differences of velocity in a numerical experiments of fluid turbulence, for typical values of r jr). 
The solid lines correspond to the Arimitsu's theory (see for details). The entire set of theoretical curves has been obtained 
with a single value qam < 1 . 

FIG. 3. Distributions of transverse momenta of hadronic jets produced in electron-positron annihilation. The solid lines 
correspond to the Bediaga-Curado-Miranda theory To each curve, a different value of q (in the range (1, 1.2)) is associated. 
The dashed line corresponds to Hagedorn's theory using BG statistics (g = 1). 



FIG. 4. Anomalous diffusion measurements of cells of Hydra viridissima The dot-dashed line corresponds to normal 
diffusion (g = 1), whereas the solid line corresponds to anomalous superdiffusion associated with q = 1.5 . 



FIG. 5. Distribution of velocities of cells of Hydra viridissima jlj]. The solid line corresponds to g = 1.5, thus coinciding 
with the value of Fig. 4. 

FIG. 6. Zipf plot (frequency of words with rank r) associated with various books as indicated on the figure (A'' is the total 
number of words; V is the vocabulary, i.e., the number of different words). See details in 

FIG. 7. Zipf plot associated with 56 books by Charles Dickens. The solid line corresponds to q = 1.9, and the crossover to 
the q = 1 regime at high rank r is visible on the figure. See details in ||l5|. 
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